Quantum and Classical Dissipative Effects on Tunnelling in Quantum Hall Bilayers 
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We discuss the interplay between transport and dissipation in quantum Hall bilayers. We show 
that quantum effects are relevant in the pseudospin picture of these systems, leading either to direct 
tunnelling currents or to quantum dissipative processes that damp oscillations around the ground 
state. These quantum effects have their origins in resonances of the classical spin system. 
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Over the past ten years, a series of elegant 
experimentsi'Si^i^ on quantum Hall bilayer systems have 
allowed investigation of the correlated many-body quan- 
tum state at total filling fraction v = 1. The bilayers con- 
sist of two closely-spaced parallel two-dimensional elec- 
tron layers in a double quantum well. If the Landau 
level fillings are i/i = = 1/2 then separate layers will 
not show the quantum Hall effect. However, Coulomb 
interactions between the layers drive a transition into a 
state in which the bilayer as a whole exhibits quantised 
Hall conductance^. Predictions^ that these states would 
undergo spontaneous breaking of a U(l) symmetry have 
been supported by experiments^. 

The ground state of the system may be viewed as 
an easy-plane ferromagnet^. There also are analogies 
with Josephson junctions^ and excitonic superfluid a^d'^ . 
The transport properties of the bilayers are rich and 
often surprising. They provide experimental evi- 
dence for the Goldstone mode associated with the bro- 
ken U(l) symmetry^ and, more recently, 'excitonic 
sup erfluidity 'i^. 

The interlayer tunnelling properties of the bi- 
layer remain a subject of considerable theoretical 
studjiii*i2iiiiiii5ii&. In a previous papeaii we introduced 
a simple thought experiment in which we investigated the 
dissipation and interlayer tunnelling properties of a clean 
bilayer system at zero temperature, and gave a brief com- 
parison of our theoretical predictions with the experimen- 
tal results. A key point is the existence of a crossoveoi^ii^ 
in behaviour at a very small bias Vq which corresponds to 
an energy of the order of the energy gap of the collective 
charge excitations in the bilayer. This gap is small for 
weak tunnelling and experiments operate principally at 
V > Vq- As mentioned in Ref.0, we predict a tunnelling 
current that varies as 1/V in that region. This negative 
differential conductance is in qualitative agreement with 
experiment. 

In this paper, we extend the discussion of Ref. ITtI 
including details of the behaviour for V < Vq, and a 
discussion of the relationship between quantum dissipa- 
tive processes and instabilities of a classical spin system. 
Section ^ contains a discussion of the model and of our 
thought experiment. In section |nl we discuss the results 
both for y > Vq, and for V < Vq- In the latter regime, 



the theoretical treatment is made rather complicated by 
the presence of an out of equilibrium Bose condensate in 
the low energy bosonic theory of the bilayer: we show 
how these issues may be resolved. In section IIIII we 
show how the quantum dissipative mechanisms in the 
system arise from instabilities of an underlying classical 
spin system. Finally we draw conclusions in section Hvl 
and identify outstanding issues in our treatment. 



I. THE MODEL 

In this section, we define our model, and describe the 
thought experiment that we will use to investigate the 
link between dissipation and tunnelling. The bilayer sys- 
tem may be cast in a pseudospin formulation? . The layer 
index of each electron defines a two dimensional Hilbert 
space: it resembles the usual spin- (1/2) degree of free- 
dom of electrons. We assume throughout that physical 
electron spin is completely polarised by the applied mag- 
netic field. In the quantum Hall state, we may work in a 
basis of localised orthogonal states in the lowest Landau 
level. Configurations in which two electrons occupy the 
same localised state represent quasiparticles, and have a 
large energy gap. Therefore, a theory of discrete spins 
(5 — 1/2) on a lattice is appropriate at low energies. 
The spins interact by an exchange interaction with an 
easy plane, and experience a field in the easy plane that 
arises from tunnelling between the layers. 

In our calculations, we generalise the model to > 
(1/2): this generahsation may be treated as a coarse- 
graining procedure. The limit of large S is also the clas- 
sical limit for the spin system. After coarse-graining, we 
replace the easy-plane exchange interaction by a combi- 
nation of isotropic exchange and an on-site anisotropy 
term. The resulting Hamiltonian is: 



H pe sr^ ^ ^ 
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where Si — Srfii — S {mf , rri^ , mf) is the spin operator on 
site z of a square lattice with spacing cq = V2ttIb where 
Ib ~ {hc/eBy/"^ is the magnetic length. The interlayer 
exchange pE and the strength of the on-site repulsion D 
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were derived from microscopic considerations by Moon 
et al& {D = Sir (31 g in the notation of that paper). The 
tunneUing between the layers enters the problem through 
AgAs: the splitting of the "bonding" and "anti-bonding" 
single-particle states in the double well. We use a Hamil- 
tonian with isotropic exchange purely for simplicity: the 
interlayer and intralayer exchange constants will differ in 
general. However, if we generalise to an exchange interac- 
tion that is anisotropic (in the spin space), then only the 
interlayer part is relevant at the large lengthscales and 
small charge imbalances that are relevant to our calcula- 
tions (we require <C 1, and psq"^ ^ D for all relevant 
wavevectors, q). Since these conditions are satisfied for 
the experimental comparisons that we make, we set both 
the interlayer and intralayer exchange constants to their 
interlayer value, pE- In our thought experiment, a gate 
is used to control the charge imbalance on the bilayer: 
this adds a term Hy — —SV^ - mf to the Hamiltonian. 

Typical values for the model parameters in physical 
bilayer systems are Ib — 20nm, Asas — 90/iK, pE — 
0.5K, D ~ 30K. 

Hamiltonians of the form |^ permit both ferromag- 
netic and quantum disordered phases at zero tempera- 
ture. At large D (or small S), the uncertainty in asso- 
ciated with ferromagnetic order means that the paramag- 
netic phase with ^ becomes energetically favoured. 
We interpret the observation of a linearly dispersing peak 
in the tunnelling conductance^ as evidence that the ex- 
perimentally accessible states have ferromagnetic order 
— the peak arises from the Goldstone mode of the sys- 
tem. In other words, we assume that the system is con- 
nected adiabatically to the large-S* limit of our model, 
ignoring the possibility of the quantum disordered phase 
that exists for D/peS^ > 1^^. 

We seek a minimal theory of the bilayer. We there- 
fore work in a clean system at T = 0. This allows us 
to calculate dissipation rates starting from purely mi- 
croscopic considerations: this is distinct from the more 
phenomenogical approach of Ref. 12 in which dissipation 
is added to the model in the form of damping terms in 
the classical equations of motion. Simulations^^ indicate 
that disorder may be relevant to the tunnelling at small 
bias voltages. While these effects are beyond the scope 
of this paper, we argue below that our framework may 
be suitable for further investigation of these effects. Pos- 
sible effects of topological defects such as merons^ are 
also neglected in our treatment, as are inhomogeneities 
in the order parameter associated with the spontaneous 
symmetry breakingi^. 

A feature of our framework is that we work on trans- 
port in an isolated bilayer. This means that our calcula- 
tions are well-controlled. Coupling the bilayer to leads, 
as in experiments, complicates the picture. We believe 
that the bulk effects that we study should dominate over 
boundary effects associated with the leads, such as those 
studied by Parkin, especially in the (experimentally rel- 
evant) geometry where the current source and drain are 
separated by a distance of the order of the sample size. 
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FIG. 1: Sketch of the constant energy trajectories on the spin 
sphere. 

Having reduced the system to a minimal model, we 
now proceed to discuss the behaviour of the spin system 
of equation 

A. Classical limit of the spin system 

In the ferromagnetic phase, we expect the qualitative 
behaviour of the system to be accessible through an ex- 
pansion about the classical limit. In this limit, an ini- 
tially uniform configuration of the spins remains uniform 
forever. Each spin precesses along one of the trajecto- 
ries of constant energy shown in figure ^ There are two 
qualitatively different types of trajectory in that figure: 
those which wind around the x-axis and those which 
wind around the z-axis. This distinction underlies the 
crossover in the tunnelling IV relation mentioned above. 
We now discuss the form of the trajectories and the clas- 
sical properties of the spin system. 

Trajectories with low energies are localised near the 
ground state, in which the spins are aligned along the x 
axis. The spins precess around that axis with a frequency 
Asw/^ where 

A,w = [Asas ( Asas + ^)]'/' (2) 

is the energy gap for spinwave excitations. These ex- 
citations describe density waves of charge imbalance 
across the bilayer. For the case Asas = these exci- 
tations form the Goldstone mode that appears when the 
symmetry around the easy plane is spontaneously bro- 
ken. In that case they disperse linearly with a veloc- 
ity V = Isi^nDpEy^^- Introducing the tunneUing Asas 
breaks the symmetry explicitly, and the gap opens up. As 
the mean magnetisation precesses around the x-axis, the 
charge imbalance on the bilayer oscillates around zero. 

For trajectories with much larger energies, the mean 
spin precesses around the maximal energy state, which 
is close to the z-axis of the spin sphere. This yields a 
Josephson-like alternating current / ~ cAsas cos{eVt/h) 
where V is the voltage across the bilayer due to capaci- 
tative charging. We stress that this is valid only for large 
charge imbalance (large V). 

A saddle point on the negative x-axis of the spin sphere 
marks the boundary between oscillations around the 
ground state and oscillations around the maximal energy 
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FIG. 2: Sketch of the charge imbalance on the bilayer (S^) in 
the thought experiment, as a function of time after the gate 
voltage is removed. 



state. The saddle point trajectory crosses the xz-plane 
at the saddle point, and at (x, z) — {S cosBq, S sinOo) 
where 6q is the angle shown in Fig. ^ It is given by 
cos 00 = 1 -2(AsAs/^). 

When the gate is used to induce a charge imbalance on 
the bilayer, the mean spin is tilted out of the easy plane 
by an angle 9, which satisfies: 



eV = Dsind + Asas tan 



(3) 



where V is the gate voltage. The saddle point then cor- 
responds to a voltage difference of 



Vo = (2A,w/e)[l + 0(AsAs/^)] 



(4) 



across the layers. (We write Vq in this form since the tun- 
nelling energy scale set by AgAS j is much smaller than the 
charging energy scale set by D.) Since the dynamics of 
this system depends strongly on whether the applied bias 
is above or below this threshold value, the calculations 
we describe below will treat these two regimes separately. 



B. A thought experiment 

We now turn to our thought experiment which will 
form the basic strategy behind our calculations. We 
imagine using a gate to induce a uniform charge imbal- 
ance on the bilayer. In the spin picture, the magnetisa- 
tion is tilted out of the easy (xy) plane with 



= sin 6* ~ eV/D 



(5) 



where V is the voltage bias applied across the bilayer and 
the second approximate equality is valid to leading order 
in (AsAs/^) and (eV/D), both of which will be small 
(recall equation ©). The magnetisation remains in the 
a;z-plane due to the tunnelling field. The bias is then 
instantaneously removed and the bilayer finds itself in a 
highly excited state. In this paper, we will calculate the 
quantum and classical dynamics of this highly excited 
system. 



From the discussion of the classical dynamics, we ex- 
pect the thought experiment to give different results ac- 
cording to whether the initial gate voltage V is above or 
below the saddle point value Vq. For small initial energy, 
V < Vq, the charge imbalance oscillates around zero: 
the effect of any dissipation will be to damp this oscil- 
lation (lower trace in figure [2Jl. Conversely, if the initial 
charge is larger, the oscillation is about a finite charge 
imbalance: dissipation leads to a decay of that charge 
imbalance, and therefore a direct current (upper trace in 
figure 121). The link between dissipation and transport is 
clear in this case, since dissipative processes lead directly 
to a direct current. 

In our calculations, we will treat these two regimes in 
separate perturbation theories. Common to both calcu- 
lations, we will expand around the classical limit using 
spinwave theory m a 1/S expansion. In a quantum the- 
ory of the spin system, the coherent oscillations of the 
charge imbalance do not remain coherent forever. This 
is different from the classical case and is consistent with 
the mixing between modes with different wavevectors ob- 
served in the exact diagonalisations of Nakajimai^. Dis- 
sipation is possible because the tunnelling term breaks 
the global spin rotation symmetry so that the long- 
wavelength modes are no longer protected from decay 
by Goldstone's theorem. We will find that the leading 
terms in the expansion about the classical limit result in 
the leading terms in the relations between the dissipation 
rate and the initial gate voltage, V. 



II. DISSIPATION RATES IN THE QUANTUM 
SYSTEM 



In this section, we calculate the dissipation rates out of 
the coherent classical oscillation of the spin system. As 
shown in figure |21 the effect of the dissipation is different, 
depending on the initial charge on the bilayer. However, 
in both cases, the dissipation rate, F is given by 



r = a 



iHt 



(6) 



where Eq measures the energy in the bosonic mode with 
momentum q and \i) is an initial state for the system. 
The procedure used to calculate the dissipation rate is 
different in each case. 

The smallest energy scale in the physical bilayer sys- 
tems is Asas j so it seems natural to treat this parameter 
perturbatively. However, the tunnelling term breaks a 
continuous symmetry of the Hamiltonian, so care must be 
taken in the limit Asas — > 0. It is clear from figure^that 
the trajectories with energies above the saddle point en- 
ergy are perturbatively connected to the AgAS = limit. 
Conversely, trajectories with very low energies are quali- 
tatively different to those in the limit of zero tunnelling: 
in this case perturbation theory in Asas is not appro- 
priate. Instead, we treat the energy of the trajectory (or 
equivalently the initial gate voltage) perturbatively. 
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We now calculate the extent to which quantum dissi- 
pative processes affect charge relaxation on the bilayer, 
or damping of charge oscillations (recall figure EJ- 



(a) 



q=0 



9 (b) 
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A. Charge relaxation: V > Vo 

We begin by considering the situation for initial ener- 
gies above the saddle point {V > Vq). This calculation 
was outlined in Ref. 17; we review the method here since 
the calculation in the following subsection is a refinement 
of the same approach. Evaluation of T requires the choice 
of an initial state, and calculation of matrix elements be- 
tween that state and final states with different values of 
Eq. We use a bosonic representation of the spin algebra, 
which defines the basis for initial and final states. We 
then identify the initial state and the dominant process 
that leads to dissipation; this allows calculation of T. 

To begin, we define our bosonic representation of the 
Hamiltonian, by using the Holstein-Primakov represen- 
tation of the spin algebra SJ + iSj = (25* — a]aj)^/^aj, 

Sj = S — ajoj. We treat tunnelling perturbatively, so 
we begin by setting Asas = and expanding around 
= {eV/D) = sin 6. The quadratic part of the Hamil- 
tonian is then easily diagonalised in the Fourier basis to 
give 



H. 



(0) _ 



{Dsint 



q#0 



(7) 



where the aq are bosonic operators ([aqiCtg/] 
describing the spinwave modes. The spinwave dispersion 
is given by 



7(9) 



4-2 cos{q^co) - 2 cos{qyCo). 



As discussed above, the tunnelling is being treated per- 
turbatively. Therefore, this dispersion relation is that of 
a model with spin rotational symmetry in the a;?;-plane, 
and Eq is gapless linear in q at long wavelengths. The 
absence of an energy gap in this approximation does not 
affect our calculation for V > Vq — 2As„/e. Observe 
that the q — mode has been singled out in the Hamil- 
tonian, and its energy is not given by the long wave- 
length limit of Eq. The quanta of this mode carry = 1 
while spin waves with finite wavevector have — 0. 
Thus the eigenspectrum of the system consists of mul- 
tiple branches: each branch has an associated S^, and 
consists of a continuum of linearly dispersing collective 
modes. We use a basis | A^, {q, q', . . .}) where N is the 
number of quanta in the q = mode (which labels the 
branch) , and the set of wavevectors indicate indicate the 
presence of collective modes with finite momenta. 

The next step is to express the initial state of the sys- 
tem in this basis. The thought experiment stipulates an 
initial state with the expectation of the spin equal to 
{S coa9,0, S sm9): this is a coherent state. However, a 



FIG. 3: Decay processes above the saddle point: (a) in the 
absence of an in-plane magnetic field, we see decay of a single 
quantum of the q — mode into a spinwave pair that satisfies 
Eq = e_q — eV/2. (b) in the presence of a in-plane field, a 
single quantum of the q — mode can decay into a spinwave 
excitation with momentum Q. 



choice of an eigenstate of with the same expectation 
of the energy leads to the same average rate. Thus we 

use\i)v>Vo = \{eV/D){SL^),{})- ^ 

The final step in the calculation is to identify the dom- 
inant process contributing to F. The dissipation arises 
from the destruction of one quantum in the q = mode, 
and the generation of multiple spinwaves during tun- 
nelling across the bilayer. To leading order in 1/S, a 
pair of spin waves is excited, each with energy eV/2 but 
with opposite momenta (Fig. . The relevant vertex in 
the interaction Hamiltonian is: 



(1) 



iSAS 



(8) 



in which Asas = Asas exp(— S'^^-^/ZJ/pi?) is a renor- 
malised tunnelling amplitude, and the vertex factor is 
given by j2q = coa9[{uq + Vq)^ + 2sin0sec^6' — {uq — 
VqY sec^ 9] where Uq and Vq are coherence factors: [uq + 
vqf = {D + PEl{q))/eq with ul-vl = 1. The 9- 
dependence arises because the x-component of the pseu- 
dospin depends on the angle 9 of the spin with the easy 
plane, as well as the azimuthal angle (j) around the plane. 
This dependence is weak for <C 1 when the charge im- 
balance is small compared to the Landau level filling. 

Hence we calculate the power dissipation F for a given 
initial voltage V, using Eq — Sqal^aq. The calculation 
proceeds according to Fermi's Golden Rule. The steady- 
state tunnelling current density at a bias V can then be 
computed from this dissipation by / = V/VL'^S. We 
find: 



L25 



V 



32TTpEllhS 



[1+xm' (9) 



where X{9) — {sec 9 — l){sec^ 9 — 2 sec 6 — 1) is small 
for a small charge imbalance {9 <C 1). The dissipation 
saturates at small initial gate voltages, leading rise to an 
increasing current as the gate voltage is reduced. 

As discussed in Ref. this form for the current is 
qualitatively consistent with the experiments of Ref. ^3, in 
that the region of negative differential conductance is ob- 
served for V > Vq ^ 6/iV. From that experimental data 
it seems that the dissipation rate does increase with in- 
creased applied voltage rather than remaining constant: 
we attribute this increase to dissipative channels not cap- 
tured within our approach. We also note that the current 
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FIG. 4: Contributions to current at finite Q 



that flows in our thought experiment is the same current 
defined in Ref. Hol but that the calculations in that paper 
neglect the 0{1/S) contribution to the current that we 
calculate. 

Quantitative comparison between theory and experi- 
ment is complicated by the fact that the area of the sam- 
ple over which tunnelling takes place remains an open 
questioriii. It may depend on the experimental geome- 
try and on spatial fluctations in the tunnelling amplitude. 
Our calculation gives a current density per unit area of 
bilayer: in order to match this data with experiments we 
must assume tunnelling taking place over an area around 
50/im^, which is signicantly smaller than the total area, 
so consistent with tunnelling taking place over localised 
regions of the bilayer. 

The calculation leading to equation may be 
straightforwardly generalised to the case in which an 
Aharonov-Bohm flux is introduced in the plane of the 
bilayer. There is a momentum scale associated with the 
new flux, given by Q = {eByd/hcY^'^ where d is the spac- 
ing between the two quantum wells, and takes a value 
around 30nm. With an appropriate gauge choice, the 
effect of the flux is to introduce a spatial dependence 
to the tunnelling term in the Hamiltonian: Asas™'^ ~* 
AsAS cos Qx + my sin Qx) . 

The result is that a new dissipative process appears in 
which a single quantum of the q = mode decays into 
a single spinwave with momentum Q. This leads to a 
contribution to the current at order 5": in other words, 
even the purely classical spin system has a decaying os- 
cillation in this case. Evaluating the leading two terms 
in the (1/(5) expansion, we arrive at 



h 



-.Q{eV - vQ){lQ) 



2peS V(el/)2 - {vQ)' 

where Q{x) is the step function. The second term is the 
two spinwave contribution that reduces to © at Q = 0. 

The (5-function in the first term indicates a resonance 
in the classical spin system when the capacitative en- 
ergy released on transferring an electron across the bi- 
layer matches the energy of the collective mode with mo- 



mentum Q. This infinitely narrow form of the peak is 
clearly unphysical. This is because we are calculating 
a transition rate to a discrete level, by specifying both 
the momentum and the energy of the final one-spinwave 
state. The final state can itself decay at higher order in 
1/5 and so the delta function will be broadened. 

Moreover, the delta function can be broadened by dis- 
order or by finite temperature. In figure^ we show the 
IV characteristic of (|10() . where an uncertainty in the 
momenta of the collective modes has been introduced to 
mimic the effect of weak disorder. Quantitative calcu- 
lation of the broadening introduced by these effects is 
beyond the scope of this paper: the degree broadening 
also controls the peak height, which again frustrates a 
detailed comparison with experiment. However, we em- 
phasised in Ref. 17 that the asymmetric lineshape of this 
peak arises from the quantum corrections that we have 
calculated: these are multi-spinwave decay channels. 



B. Damping of oscillations: V < Vo 

We now turn to initial charging voltages V < Vq, in 
which the effect of dissipation is to damp oscillations in 
the charge imbalance on the bilayer. This behaviour was 
shown in the lower trace of figure |21 We note that this 
case is different from y > Vq in that our thought ex- 
periment does not result in a state with a steady direct 
current. For that reason we cannot directly address the 

IV relation at these small biases: the strong feature in 
the differential conductance observed experimentally^ at 
small V is beyond the scope of this paper. However, the 
coherent oscillations of charge between the layers that are 
predicted by the classical limit of the spin system are not 
observed in experiments. The intrinsic dissipation rate 
that we calculate sets a minimal decay rate for these os- 
cillations. Adding disorder or inelastic scattering to the 
model might increase this decay rate: we consider only 
the clean system in this paper, but the framework that 
we now describe could be extended to consider this case. 

The calculation of T is slightly more complex than for 

V <Vo, but proceeds along the same lines. The result is 
a slow (power law) decay of the oscillations predicted by 
the classical limit of the pseudospin system. 

In order to obtain a bosonic theory, we again use 
the Holstein-Primakov representation. However, we may 
choose either the x axis or z axis as the principal axis 
for our bosonic theory. Choosing the x axis is natural 
since that is the direction of the mean magnetisation. 
On the other hand, choosing the z axis has the advan- 
tage that no interactions arise from the anisotropy term 
in the Hamiltonian. We have checked that the results we 
present below are independent of this choice. 

Let us consider first the spinwave spectrum cuq. The 
non-interacting part of the Hamiltonian for the spinwaves 
takes the form 



(11) 
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where /3| and (3q are the creation and annihilation oper- 
ators for a spinwave with momentum q and 

LOq = [(AsAS + pEl{q)){D + AsAS + pEliq))]^'"^ ■ (12) 

where 7(9) was defined in section III Al above. Since we 
are treating the tunnelling non-perturbatively, the spin- 
waves now have an energy gap of Ag^ — [Asas(^ + 

ASAS)]^/^ 

The non-interacting eigenspectrum has multiple 
branches, labelled by the total number of spinwave 
quanta. This should be contrasted with the high-bias 
regime discussed in the previous section where there is 
a qualitative difference between the g = and finite-q 
modes and the branches are labelled by the number of 
q — Q quanta. However, for consistency, we continue to 
use the notation in which | A'o, {q, . . .}) describes a state 
with iVp quanta in the q = Q mode, as well as a set of 
collective excitations with finite momenta. (We omit the 
set of finite momenta if it is empty.) 

To set up the calculation for our thought experiment, 
the initial state of the system is chosen to be the ground 
state of the bilayer in the presence of the gate voltage 
V . In terms of the creation operator /3^^q for the q — Q 
spinwave, this is approximated by a coherent state, which 
can be written as 



\vv<v„ 



A^exp 



|0) 



(13) 



where Wq — ■\/wq/(AsAS + PEjiq)) is a coherence factor, 
and JV = cxp[—9^w^^qSL'^/4:Cq] enforces normalisation. 
The expectation of the occupancy of the q — mode in 
this state is 



iVo = noL^ = iSLy2cl)0W{D/AsAs) 



1 



(14) 



We observe that the occupancy of the bosonic mode with 
q = is proportional both to the area of the system and 
to S: the bosons in this mode are Bose-condensed with 
areal density no. We will see below that interactions 
within the condensate must be treated carefully. 

We can now discuss the Hamiltonian for our calcula- 
tion in this low-bias regime. Our perturbative parameter 
for y < Vb is {uq/S). It is independent of S and does 
not vanish in the classical limit. We will calculate the 
dissipation rate F to leading order in this parameter. We 
give the full details of the Hamiltonian in the Appendix. 
It can be written in the form 



diss 



(15) 



H[^^ describes the interactions within the condensate: 



L^S 



74<i/3(l/3o + 74e/3o + h.c. 



(16) 



and ffdiss contains dissipative processes involving finite- 
momentum spinwaves 

Hdiss — 



+lic{P% + /3lg/3-q)/3o/?o + h.c. 



(17) 



where we have included only terms that affect the dissi- 
pation rate, F, to leading order in [1/ S). The vertex fac- 
tors, 7's, are functions of q which depend on the choice of 
principal axis in the Holstein-Primakoff transformation. 
They are independent of L and are finite in the large-S' 
limit. 

Before discussing the dissipative processes, we note 
that we have isolated the interactions within the conden- 
sate in iJ^j because these processes do not directly con- 
tribute to the dissipation of energy from the q = modes 
into finite- momentum spinwaves. However, at first sight, 
they do appear to contribute to the decay of our initial 
state. Moreover, the rate for the process is apparently 
divergent in both the thermodynamic (large L) and clas- 
sical (large S) limits! For instance, for the process 74e, 
Fermi's Golden Rule would give a transition amplitude 
of 



(18) 



This seems to suggest that the condensate with A^o <Z = 
modes decays immediately without the generation of 
finite-q modes. This divergence of the Fermi's Golden 
Rule amplitude originates from the macroscopic occupa- 
tion of the q = Q mode. 

In fact, these condensate interactions do not cause de- 
cay within the condensate. Instead, they contribute to a 
coherent non-linear evolution of all the condensed q — Q 
bosons: this effect is beyond the perturbative framework 
of Fermi's Golden Rule. We will see in the next section 
that this evolution is part of the semiclassical spin dy- 
namics. To avoid divergences, we need to modify the 
interaction-picture operators to take into account the 
condensate interactions: 0(t) = exp[ii/i]0 exp[— zHt] 
where H = Hq -t- Hf^^^ contains the non-quadratic inter- 
action terms in [3q=a. This removes the divergent terms 
from the Fermi's Golden Rule calculation, but we pay the 
price of a non-trivial time dependence for the interaction- 
picture operators. For details of the modifications, in- 
cluding some comments on its physical interpretation, 
see appendix IXI 

The effect of the non-linear evolution of the q = Q mode 
is illustrated by the annihilation operator /3q=o (^) for the 
g = mode. For example, we have: 



(A^o-3|/3,=o(i)|A^o) = 



LP'S 2Asw 



{N^ - 3|/3, 



Thus Pq={)[t) has a finite amplitude for destroying more 
than one particle. (There is an analogous amplitude 
for particle creation as well.) As a result of this non- 
trivial time dependence, the 74a term in -ffdiss may now 
contribute to the decay of the initial state, giving rise 
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q=0- 
q=0- 
q=0- 
q=0- 



^ q 
^ -q 



FIG. 5: Schematic decay process for V < Vo- Energy conser- 
vation requires 20;, — 4Aaw 



to a contribution to the dissipation rate proportional to 

74a74rf- 

We now have all the ingredients required to calculate 
the dissipation rate. We note that the dissipative part, 
Hdiss, contains only terms of even order in the bosonic op- 
erators. Together with energy-momentum conservation, 
we can see that the simplest transition in which energy 
is transferred out of the q = mode involves the anni- 
hilation of four quanta of that mode, combined with the 
creation of a pair of spinwaves with equal and opposite 
momenta: 



|iVo 



\No-A,{q,~q}) 



(20) 



where the momentum q has to satisfy energy conserva- 
tion: 4Asw = 2eq. This kinematic constraint leads to an 
intrinsic dissipation which depends on the fourth power 
of the density of spinwaves. At low bias, this density is 
small and so the intrinsic dissipation is rather weak. 

As expected from kinematics, the dominant processes 
contributing to dissipation are of the form shown in fig- 
ure|3| The modified interaction picture yields an effective 
matrix element for the transition, given by 

7/,; = A;;„^(noCo/5)^(74f,74e + 74a74c " 74a74d " Asw76) 

(21) 

where the coefhcients (74a, 74&, ■ • ■) were defined in equa- 
tion (|17|l as the vertex factors for various processes in- 
volving spin wave quanta. The different terms in 7/^ cor- 
respond to different coherent routes for destroying 4 g = 
spin wave quanta and generating a pair of finite-g spin 
waves. This can be achieved with a direct vertex (76), or 
through an intermediate virtual state jA^o ^ 2, {q, — q}) 
(74q followed by 74c). The 74^ term in i/diss also con- 
tributes to the decay through the non-linear component 
of Pq^Q from ((injl giving the term 74a74rf- The 74fc74e 
term arises in a similar way. 

Substituting the relevant coefficients into (|21|) , and de- 
noting the (small) number (Asas/D) by x, the result is 



7/, = A^nocJS) 128(1 + :.)3 

(22) 

We have checked that this quantity is independent of 
the choice of quantisation axis. The density of states for 
spin wave pairs, in which the spin waves have equal and 
opposite momentum, and each has energy e is 
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The dissipation rate is now simply given by F = 
87rL2Asw7/iff(2Asw)/?i, and we arrive at our result for 
the dissipation rate in the low-bias regime {V < Vq): 



V 



^"^^^:[l+0(AsAs/i?)] (24) 



5127r/|/9£;i:>2 I 



We can express this as the rate at which the number 
of g = modes is decreasing: dNo/dt = — r(V)/4Asw 
However, Aq is related to the instantaneous bias V by 
(|14|l . The instantaneous bias is in turn related to the z- 
magnetisation by lO): V oa niz ^ 6 dX low bias. In terms 
of the tilt angle, we see that dd^/dt cx ~9^. Hence the 
amplitude. A, of the oscillating charge imbalance on the 
bilayer decays slowly in time, according to 



A{t) - 1/^1/6 



(25) 



47rp£;C^ y/D^ + 4£2 



(23) 



This is weak dissipation, but we believe that it is the 
first microscopic calculation of an intrinsic dissipation 
rate in these systems. The existence of such mechanisms 
is an important issue because coherent oscillations of the 
charge are not observed in physical bilayer systems. In- 
troducing disorder or inelastic scattering from topologi- 
cal excitations might relax the kinematic constraint that 
lead to the weak dissipation. SimulationE^"* indicate that 
these effects are important at small energies: the frame- 
work of the modified interaction picture may be a useful 
framework in which to investigate these effects. 

Equations ((Tni) and represent the main re- 
sults of this section. We have calculated dissipation 
rates in the three regimes of our thought experiment that 
are most relevant to experimental systems. As we com- 
mented in Ref.0, the results agree qualitatively with ex- 
periment; more quantitative agreement may require the 
taking into account both disorder, and the effects of finite 
temperature. 



III. CLASSICAL ORIGIN OF THE 
DISSIPATION 

In this section we show that the dissipation mecha- 
nisms in the bilayer have their origins in dynamical in- 
stabilities of the classical spin system with a Hamilto- 
nian given by (Q). We consider the case < Vq, and 
emphasise the strong links between the calculation of 
the strength of the classical resonances and the quan- 
tum calculation above. We emphasised in the previous 
section that an initially coherent oscillation remains co- 
herent forever in the classical system. However, we show 
here that there are instabilities of the classical dynamics 
whereby an initially very small spatial modulation of the 
spin field will grow resonantly, leading to dissipation out 
of the coherent oscillation. In the quantum system, zero 
point fiuctuations are sufficient to seed this instability: 
this is the cause of the dissipation rates calculated in the 
previous section. 
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The classical equations of motion of the spin system 
are described by the Landau-Lifshitz equation: 



(26) 



abc 



where Si = {Sl,Sf,Sf) is the spin on site i and e' 
is the completely antisymmetric tensor. We parametrise 
the magnetisation in terms of a mean value and a weak 



spatial modulation : rfij — mg + (mqC 



c.c). As- 



suming that the spatial modulation is small compared to 
the mean value, we may linearise the equation of motion 
for mo, arriving at 



{hs) 

(hS) 







dt 
drfiq 

"df 



Too X [AsAsea; - De^TO^] 



(27) 



TOq X [AsAsei - ^'e^TOol 

+TOo X [-DelTOq - pE-/iq)mq] (28) 



where ex,y,z are the unit vectors in the a;,?;,z-directions 
in spin space. 

Let us consider first the time evolution of the uni- 
form mode TOo(t). Suppose we have an initial state with 
a small component in the S'^-direction: TOg(i = 0) = 
{y/l — Zq, 0, zq) with zq <^ 1. We will use zq as our per- 
turbative parameter. (Compare with the quantum calcu- 
lation of section III Bl where the perturbative parameter 
was the energy in the q = mode.) 

At very small amplitudes, the driving mode is har- 
monic: the spin precesses along elliptical trajectories 
with frequency Agw/^- The perturbation expansion over 
Zq is straightforward, and is equivalent to solving for the 
motion of the operator f3q=o(t) in the modified interac- 
tion picture (see appendix IXIl . The time dependence of 
the z-component of the uniform mode ttiq is 



TOg(t) = zocos{Aswt/h) 



D4 
64A2 . 



cos(3A,wt/?i) - cos(A,wV?i) + 0(^0X29) 



where A^w = Asw[l - {Dza/A^A^^f + 0{z^)\. The other 
components of tuq are easily obtained from Asas^o ~ 
{hS){dm'^/dt) and {m%f = 1 - (to^)^ - {mlf . 

Let us now turn to the spatial modulations rhq{t). It 
is simple to verify that |too| and (too • to^) are constants 
of the motion. To make this explicit, we write 



mq =P2q(™o X el) +Piq{rno x (too x el)) 



(30) 



This is the equivalent of casting the spin system (with 
three components per site) in terms of bosons (one com- 
plex or two real degrees of freedom) via the Holstcin- 
Primakov representation. In this representation, the lin- 
earised equation of motion for the spatial modulation is: 



i^S)- 



PlQ 



dt \ P2q 



1 - 



Mdq Maq 
~Mbq Mdq 



Plq 
P2q 



(31) 



2A-e 




E-2A 



FIG. 6: 'Energy level' picture of the classical calculation. Res- 
onance occurs when initial (left) and final (right) states be- 
come degenerate. 



AsASTOi5 + 

\- PEiiq) + 



where Mdq = AsAsmgTOg/(l- (TOg)^), M^q = 
PEliq) and Mbq = [AsAsmg/(l - {m'of)] 

Z?(l-(TOg)2). 

Consider first the limit of small amplitude zq ~* 0- 
Then, Md = 0, Ma = Asas+ pEl{q) and Mb = PEl{q) + 
D. We can see that each g-mode is harmonic with a 
natural frequency ujq/h = (MaqMdq)^^'^ /h which is the 
same as the dispersion relation given in H12|l . These are 
the classical spinwaves of the system. 

However, at larger amplitudes zq, we cannot ignore the 
driving of these spinwaves by the uniform mode. This will 
provide the energy to amplify any small spatial modula- 
tions. Substituting mo{t) into |(2U, we can see that these 
matrix elements are periodic in time. They take the form 

Mq''{t) = ^^tn e^^/'')"'^™* (32) 

n 

We may therefore find solutions of the Bloch form 



Plq 
P2q 



ik{q)t ( Ulq{t) 
U2q{t) 



(33) 



where ui and U2 are periodic in t with period (27r?i/ Agw)- 
The relation between the wavevector of the spatial mod- 
ulation, q, and the Bloch wavenumber, k, is a complex 
band structure, in which complex wavenumbers represent 
resonantly growing (or exponentially decaying) solutions. 
The relation k[q) is obtained by solving, for any given 
spatial mode q, an equation of the form 



det 

ab.mn 



{[/c(q) + 2nAsw]<5„ 



ab 



iM. 



ab 



= (34) 



where the determinant is of an (infinite dimensional) ma- 
trix operating in the product space of Fourier components 
(mn) and components {ab) of the vector p. 

This calculation closely resembles the calculation of the 
quantum matrix element between initial and final states: 
see figure El The non-interacting {zq = 0) problem de- 
fines energy levels which we may identify as initial and 
final states. At finite zq, we treat the interactions be- 
tween the states perturbatively: the resonance condition 
occurs when initial and final states become close in en- 
ergy (as in the quantum case). 

The resulting complex band structure for spatial 
wavevectors around ujq = 2Asw is shown in figure \7\ 
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FIG. 7: Form of the complex band structure at the resonance 
{ujq = 2Asw). We abbreviate the driving frequency Asw to A. 

There is a resonance at 2ujq = 4Asw At the centre of 
the resonance, the imaginary part of the Bloch wavenum- 
ber, A, coincides exactly with the matrix element for the 
quantum decay process. To be precise, 

A = (A,w/?i)(^o/^)'(Asw/2AsAs)' x 

3 -Ax- 2xVl + I6x + 16x2 - Sx^ 

128(1 + ^^^^ 

where we have again abbreviated (Asas/^) x- Identi- 
fying with (zq/S), and using equation H14|l in conjunc- 
tion with equation 122|) . we see that 

A = Jf^/h (36) 

Equation is strong evidence that the quantum de- 
cay process is linked to the instability of the classical 
system at the same wavevector. (Note that even the non- 
trivial dependence on As as / D contained in the factor jfi 
is identical in both cases). 

To make this correspondence closer, we demonstrate 
that with an appropriate choice of boundary conditions, 
the rate at which the classical oscillation decays is ap- 
proximately equal to the rate given in (|24() . An outline 
of this calculation is given in appendix^ The boundary 
conditions of the classical calculation have been chosen 
so as to mimic the quantum system. To set the bound- 
ary conditions, the initial energy in the driven mode is 
of the order of zero point energy in that mode, and the 
initial phase difference between driven and driving mode 
is chosen such that the energy can never drop below its 
zero point value. 

The significance of the result (|B9|I is that the dissipa- 
tion rate in the classical system takes the universal form: 

Fci - 4^i;(0).9c(2A,w)A2C(At) (37) 

where E{0) is the initial energy in the spinwave modes 
with frequency (2Asw/?i), .9c is the classical density of 
states at that frequency, and the form oiC(t) is shown in 
figure 

Figure IHl shows the similarities and differences between 
the classical and quantum dissipation rates. For the clas- 
sical trace, we used E{0) = 2Asw, which is of the order 



FIG. 8: Classical (solid) and quantum (dotted) dissipation 
rates, in units of the (constant) quantum rate. The time 
dependence of the classical rate follows the universal function 
C{t): see (I37II and appendix f51 

of the zero point energy in the relevant spinwave mode. 
At short times {Xt ^1), the classical system has a rapid 
transfer of energy into modes far away from the reso- 
nance. Once this transient behaviour is over, the quan- 
tum and classical rates are both approximately constant 
in time at intermediate times. The two results diverge 
again for long times At 3> 1. However, the quantum 
calculation is non-degenerate perturbation theory, and is 
therefore valid only up to times Xt ^ 1 and so we do not 
expect the agreement to continue beyond this time scale. 
We therefore emphasise that the form of the dissipation 
rate calculated in appendix IbI and plotted in figure |H1 is 
a universal result for dissipation out of a periodic oscil- 
lation into a bath of "driven" oscillators (as long as the 
equations of motion of the driven oscillators may be lin- 
earised). This may be regarded as the classical analogue 
of Fermi's Golden Rule. 



IV. CONCLUSION 

We have used a thought experiment on an isolated, 
clean system at zero temperature to investigate the link 
between tunnelling transport and dissipation in quantum 
Hall bilayers. The thought experiment allows us to iden- 
tify a crossover in the behaviour of the system that may 
be underlying the small bias crossover in the tunnelling 
IV characteristic. Our theory predicts an IV character- 
istic that is in qualitative agreement with experiment at 
large bias. 

In the very low bias regime we have identified an in- 
trinsic source of dissipation in the system, and showed 
that its origin is in a dynamic instability of a classical 
spin system. 

In the light of these results, two main questions remain. 
Firstly, we were unable to address the IV characteristic 
at small biases, which is the position of a prominent fea- 
ture in the conductance. An extension to the method 
so that a steady-state current can be forced through the 
bilayer would be required if this issue is to be resolved. 
Secondly, the effects of disorder on the model have not 
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been considered, although this thought experiment and 
theoretical framework does seem appropriate for investi- 
gating that problem. 
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APPENDIX A: 



MODIFIED INTERACTION 
PICTURE 



In this appendix, we show how modifying the usual 
interaction picture allows us to avoid divergences associ- 
ated with interactions between quanta of the macroscopi- 
cally occupied q = mode in the low-bias regime (section 
IIIBII . We begin by describing our modifications, before 
commenting briefly on their physical interpretation. 

For our calculation of the dissipation rate ^ for the 
low-bias regime in section FlI Bl we wish to evaluate ma- 
trix elements of the form 



X-(0,{g,-g}|A 



N-m^iHt 



Eqe 



-iHt(aUN 



{Pir\o,{}) (Ai) 



where the operator Eq measures the energy in the spin- 
wave mode with finite momentum q and H is the Hamil- 
tonian from (|15|) . is a macroscopically large number 
and m is a number of order unity. For convenience, we 
use the notation /3o = Pq=o for the annihilation operator 
of the q = mode. 

In the usual interaction picture, we would use time- 
dependent operators 0{t) — e*^°*Oe~*^"* and the 
wavefunctions have a time dependence of \ip{t)) ~ 
IV'(O)). We would write the matrix element 



aiHat 



iHt 



(0, {q, -g}|/3^'^"e*"*e-^^°'£;,e^««*e-^^*(/3T)'" |0, {}) 

(A2) 

where is the quadratic Hamiltonian (|lll) that com- 
mutes with Eq. The product lJ{t) — ^iHot^-iHt ^^^^ 
be written as a time-ordered exponential: U{t) =^ 
rexp[— i At' Hi{t'y\ which has a perturbative expan- 
sion in terms of time integrals of the operator Hi (t) = 

Our modification to the interaction picture replaces the 
time evolution e^J-^°*' in the definition of the interaction 
picture with e~'^* where H includes the non-quadratic 
terms from H^^^ (I16|l . but maintaining the commutation 
relation H] — for all modes except q — Q. Thus H 
includes all terms that are functions only of (/3o,/3o): we 
have 



H = 



E 

1 



[74d(/3^/3o' + h.c.) + 74e(/3o^ + h.c.)](A3) 



where the omitted terms have coefficients that are 
0[{L^S)~^]. The operator Hi{t) must now be evaluated 
in a perturbative series before the expansion of the expo- 
nential in U(t). This procedure requires the (non-trivial) 
time dependence of the operator /3o which is given (for 
TO ^ 0) by 

{N-m\e'"'l3oe-'"*\N) = ^ 5,n-ie-'^'-' 
{N-m\[P^,Ht^,]\N) 



'NO{N/L^Sf (A4) 

These matrix elements converge even in the classical (or 
thermodynamic) limit, in which N ^ oo but NS^^ and 
NL~'^ are constants. 

The purpose of this calculation is to evaluate a decay 
rate into a continuum of final states. To interpret the 
modified interaction picture in this light, observe that 
we can write the matrix element, M in the form 

(A5) 



M 



Y,{f\Eq\f)\{f\e- 
\f) 



where \i) and |/) arc initial and final states. Then, writ- 



ing I/) 
that 



JHt 



| g), the fact that H and Eq commute means 



M 



Is) 



{g\Eq\g)\{g\U{tm' 



(A6) 



Note that the state \g) — e^*^'|/) is obtained from |/) 
by evolving the state |/) in parallel with the initial state, 
taking into account the evolution of the condensate itself, 
but ignoring interactions that involve modes with finite 
wavevectors. This does not affect our results because 
the quantity of interest, Eq, is independent of the occu- 
pancy of the condensate: formally {g\Eq\g) = {f\Eq\f). 
In other words, this is a unitary change of basis for our 
subspace of final states of the form |7V, {q, q', ...}). 

The advantage of using the \g) basis is that it allows us 
to use conventional perturbation theory for the evolution 
operator U{t). Otherwise, we find divergences when eval- 
uating matrix elements to individual final states in the 
I/) basis. Physically, these divergences arise from the 
fact that the particles in the condensate have evolved 
in time and so the initial condensate |A^o) has become 
strongly admixed with states with different numbers of 
(7 = modes: |A^o + q})- The \g) basis tracks 

this evolution and is therefore the natural basis for the 
calculation. 

Equations lASp and (|A6|) define a convergent pertur- 
bative expansion for the dissipation rate F, defined in l^: 
this leads to equation (|24|l . 



APPENDIX B: DISSIPATION RATES IN 
CLASSICAL SYSTEMS 

In this appendix we discuss the extent to which a 
broad class of quantum problems with solutions based 
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on Fermi's Golden Rule may be addressed within a clas- 
sical framework. We investigate the decay of a coherent 
oscillation by dissipation into an environment modelled 
by a 'bath' of harmonic oscillators. The only quantum 
ingredient will be the zero-point fluctuations of this bath. 

To illustrate the point we use a simple model system 
with a Hamiltonian of the form 



H 



1 



[pin + ^l^l) + AXn^xly^' 



(Bl) 

where the coordinate y describes the coherent (driving) 
oscillation and the coordinates x„ describe the bath of 
driven oscillators. We work perturbatively around A = 
which is the limit of independent oscillators in which 
the coherent oscillation remains coherent for ever. The 
energy in the driven oscillator is initially much greater 
than all other energy scales, so y ~ cos fit, and the Xn 
obey Mathieu's equation^'' for small a;„. 



{dt +ujI + X sin 217t)a;„ ~ 



(B2) 



For each n, solutions can be written in the Bloch form 
Xn{t) = e^''*Un{t) where u„(<) has a period of 27r/ri. 

The interesting behaviour occurs when (w„/r2) is very 
close to an integer. The Bloch wavenumber k becomes 
complex and gaps open up. Near w„ = fi, it can be 
shown that 



(fc -I- nf ^ K - nf - xl 



(B3) 



Thus there is a gap at the edge of the Brillouin zone: 
Re(fc) = — f2. For |a;„ — i7| > A„, fc is real and the 
solution for a given n-mode is 



where R and </> are determined by boundary conditions. 
For \uJn — r2| < A„, k is complex. Writing k = i{k — VL) 
(which is real) and 20 — arg(iK — LOn -f fi), the solution is 



x^ = Ae"* cos(m + 9)- Be^"* cosint - 



(B5) 



where A and B are determined by initial conditions. 

In our quantum calculation based on Fermi's Golden 
Rule, we calculated the energy transferred to the driven 
oscillators, using for the energy of each mode the energy 
given by the non-interaction Hamiltonian. The analogous 
quantity here would be the dissipation rate 



T,,^dtY.{l/2){pl^+u;lxl). 



(B6) 



Using the solutions above, the energy in each driven 
mode with real k is 



A„ coB[2{k + n)t + 2(j)'] 

Wn — Xn COS 2(j)' 



(B7) 



with w„ = (A^j + \k + r2p)i/2^ The angle 4>' is equal to 4> 
for > fi: if w < ri then (j)' — 0-1- (7r/2). Similarly, if k 
is complex then we have 



A„(cosh2Kt + sinh2Kt sin0") — Wn cost 



A. 



= E{Q) 

(B8) 

where the phase angle (p" depends on the ratio of A and 
B in equation (|B5p . 

For correspondence between the classical and quan- 
tum calculations, the initial energy in the driven oscil- 
lators should be of the order of the zero point value of 
the quantum calculation. Further, we choose the relative 
phase between driven and driving oscillators such that 
the energy in the classical oscillators can never be re- 
duced below the same zero point value (since that would 
violate the uncertainty principle). 

Finally the result is that 



E{0)X^g^{n)C{Xt) 



(B9) 



where gd'^) is the classical density of states : ^ 
J duj gc{u}) in the limit in the frequencies in the oscillator 
bath form a continumir. The universal function C{Xt) 
was plotted in figure |S1 It takes the form 



C{x) 



{1/2ttx) - 
l + 0{x) 



(1/4) x < 1 
a; > 1 



(BIO) 



The diverging rate at small x is cut off by the finite width 
of the oscillator bath. At small times there is a rapid 
transfer of a small amount of energy into modes far from 
the resonance. 

The rate calculated in (|B9|I is compared with the dissi- 
pation rate for a quantum calculation on the same system 
in figure |S1 The agreement breaks down at long times 
since the classical calculation corresponds to degenerate 
perturbation theory for states near the resonance, and is 
able to capture the resonant growth of the energy. At 
small times there is a transient effect that is unique to 
the classical system. However, at intermediate times then 
choosing Eq = fiUn gives good agreement between quan- 
tum and classical rates. This is in accordance with the 
correspondence principle, as expected for systems with 
linear equations of motion. 

The applicability of this treatment to more compli- 
cated models than that of equation ljBl|) is not obvious 
at first sight. However, the forms for the energy in the 
driven modes, HB7|I and HB8(I . are in fact universalSi for 
problems in which the equation of motion of the driven 
modes is linearisable (that is, problems in which the pe- 
riodic forcing multiplies a coordinate as in ljB2|l . so that 
the forcing vanishes as the energy in the driven oscillator 
gets small). In this general problem the parameter A in 
(|B9|) is equal to the maximal value of the imaginary part 
of the Bloch wavenumber. Neglecting the effect of the 
initial transient on the upper cutoff, ljB9|l is then a uni- 
versal function of Xt, so the dissipation rate associated 
with resonance is completely characterised by the A and 
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the density of driven oscillators gd^)- This is the gen- and shows that the quantum decay channels have corre- 
eral result applicable to the spin system of equation sponding instabilities in the classical dynamics. 
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